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Abstract 

Bog  Spicebush  (Lindera  subcoriacea)  is  a  U.S.  Department  of  Defense 
(DoD)  Species  at  Risk  (SAR)  that  has  recently  been  proposed  for  listing 
under  the  Endangered  Species  Act.  Roughly  60%  of  all  known  Bog  Spice¬ 
bush  populations  are  found  on  five  DoD  installations:  Fort  Bragg,  Camp 
Mackall,  Fort  J  ackson,  Eglin  Air  Force  Base  (AFB)  and  Camp  Shelby  J  oint 
Forces  Training  Center  (J  FTC) .  This  recently  described  species  has  been 
reported  to  occupy  a  variety  of  plant  communities  that  have  varying  dis¬ 
turbance  dependence,  suggesting  the  habitat  suitability  of  occupied  sites  is 
spatially  and  temporally  dynamic.  Additional  information  about  Bog 
Spicebush  habitat  reguirements  and  potential  distribution  is  needed  to  in¬ 
form  assessments  of  the  species'  conservation  status  and  management 
needs.  This  project  used  MaxEnt,  a  widely  applied  species  distribution 
modeling  approach  to  generate  a  range- wide  habitat  suitability  map.  This 
map  was  also  used  to  identify:  ( 1)  sites  warranting  targeted  surveys  for 
novel  populations,  (2)  suitable  habitat  for  conservation  and  management, 
and  (3)  sites  suitable  for  population  (re)introduction. 
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1  Introduction 

1.1  Background 

Land  management  agencies  increasingly  recognize  the  utility  of  species  dis¬ 
tribution  modeling  (SDM)  for  rare  plant  conservation  planning  and  man¬ 
agement  (e.g.,  Elith  and  Burgman  2002,  Williams  et  al.  2009,  Gogol- Proku- 
rat  2011).  Knowledge  about  species'  habitat  requirements  and  the  spatial 
distribution  of  available  habitat  are  critical  for  designing  and  implementing 
many  conservation  actions,  but  oftentimes  this  knowledge  is  incomplete  for 
rare  plants.  Species  distribution  modeling  can  provide  insights  about  habi¬ 
tat  requirements  and  availability,  which  in  turn  can  be  used  to: 

L  Guide  targeted  surveys  that  enhance  discovery  of  previously  unknown 
populations 

2.  Prioritize  land  management  and  acquisition  initiatives  to  better  support 
species  conservation  and  recovery 

3.  Foster  additional  opportunities  for  sharing  conservation  responsibilities 
among  Federal,  state,  and  Nongovernmental  Organization  (NGO)  land 
management  partners 

4.  Confidently  select  sites  suitable  for  population  (re)introductions. 

Bog  Spicebush  (Lindera  subcoriacea)  is  a  rare,  dioecious,  avian  dispersed, 
woody  shrub  of  southeastern  wetlands.  The  species  was  only  identified 
and  described  in  1983,  based  on  leaf  morphology  and  odor;  previously  it 
had  been  included  under  the  more  widespread  Northern  Spicebush  (Lin¬ 
dera  benzoin)  (Wofford  1983).  The  known  range  of  Bog  Spicebush  spans 
seven  southeastern  states,  with  a  large  number  of  populations  occurring 
only  in  North  Carolina  and  Mississippi.  It  is  reported  to  occupy  habitats 
having  varying  disturbance  dependence,  namely  Atlantic  white  cedar  for¬ 
ests,  coastal  plain  small  stream  swamps,  streamhead  pocosins,  and 
sandhill  seeps  in  North  and  South  Carolina  (Some,  Gray,  and  Crutchfield 
2006;  Wall  et  al.  2013),  and  the  wettest  portions  of  sphagnous  bogs  in 
Mississippi  and  Florida  (Gordon,  J  ones,  and  Wiseman  1986) . 

Published  studies  for  Northern  Spicebush  (Qpollini,  Wallace- Senft,  and 
Whigham  1994;  McEuen  and  Curran  2006)  and  preliminary  demographic 
work  (Wall  and  Hohmann,  unpublished  data  2016)  suggest  that  Bog  Spice- 
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bush  populations  are  likely  dynamic,  with  metapopulation  processes  affect¬ 
ing  recruitment  into  sites  that  vary  in  suitability  over  time  and  space.  Alt¬ 
hough  Bog  Spicebush  has  roughly  the  same  state  and  global  conservation 
ranks  as  the  Federally  endangered  Pondberry  (L.  melissifolia)  (NatureServe 
20 13),  comparatively  limited  information  is  available  in  the  published  liter¬ 
ature  regarding  the  species'  biology,  population  dynamics,  habitat  reguire- 
ments,  or  distribution. 

Bog  Spicebush  was  included  in  a  recent  court  settlement  that  mandates  the 
U.S.  Fish  and  Wildlife  Service  (USFWS)  to  make  a  determination  of  listing 
under  the  Endangered  Species  Act  (ESA) .  In  addition,  it  is  identified  by  the 
U.S.  Department  of  Defense  (DoD)  as  a  Species  at  Risk  (SAR)  (NatureServe 
20 11, 20 14) .  The  majority  ( -60%)  of  all  known  Bog  Spicebush  populations 
occur  on  five  DoD  installations,  with  roughly  47%  and  10%  on  Fort 
Bragg/  Camp  Mackall  (NC)  and  Camp  ShelbyJ  FTC  (MS),  respectively.  The 
species  is  also  found  on  Fort  J  ackson  (SC)  and  Eglin  Air  Force  Base  (AFB) 
(FL).  Unfortunately,  long-term  prospects  for  Bog  Spicebush  are  not  opti¬ 
mistic  without  concerted  management  intervention.  It  is  projected  that  the 
number  of  populations  in  North  Carolina,  which  were  recently  estimated  to 
have  suffered  a  28%  decline  over  the  last  30  years,  will  undergo  additional 
declines  given  that  many  of  the  small  populations  (20%  had  <2  individuals) 
are  at  risk  of  extirpation  (Wall  et  al.  20 13) . 

The  rarity,  small  size,  and  declining  number  of  Bog  Spicebush  populations, 
as  well  as  limited  information  about  the  species'  ecology,  indicate  species 
distribution  modeling  would  be  invaluable  for  informing  the  listing  deci¬ 
sion,  conservation  planning,  and  proactive  management. 

1.2  Objectives 

The  objectives  of  this  effort  were  to: 

1  Model  the  potential  distribution  of  the  species 

2.  Provide  land  managers  with  a  three-part  strategy  for  Bog  Spicebush  con¬ 
servation  by: 

a.  identifying  sites  warranting  additional  survey  effort 

b.  prioritizing  sites  for  protection/ management 

c.  identifying  sites  suitable  for  potential  (re)introduction  efforts. 
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1.3  Approach 

The  objectives  of  this  work  were  accomplished  in  three  primary  tasks: 

(1)  data  acquisition  and  preparation,  (2)  model  development,  and  (3)  con¬ 
servation  assessments,  based  on  model  outputs. 

Secondary  tasks  associated  with  model  development  included:  (1)  generat¬ 
ing  training  and  testing  datasets,  (2)  defining  the  spatial  extent,  spatial 
thinning  distance  and  number  of  pseudoabsences,  (3)  reducing  the  num¬ 
ber  of  environmental  predictors  and  assessing  overfitting,  and  (4)  running 
the  finalized  model. 

1.4  Scope 

This  effort  modeled  Lindera  subcoriacea  habitat  suitability  across  seven 
southeastern  states  using  occurrence  data  and  a  diverse  suite  of  environ¬ 
mental  predictors,  including:  climate,  topography,  fire  regime,  vegetation, 
and  soils.  Conservation  assessments  based  on  the  model  output  empha¬ 
sized  not  only  DoD  installations  known  to  support  populations  of  L.  sub¬ 
coriacea,  but  also  the  myriad  public  and  private  land  managers  that  might 
serve  as  potential  conservation  partners. 
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2  Methods 

2.1  Model  Development 

2.1.1  Species  occurrence  data 

Species  occurrence  data  were  acquired  as  point  and  polygon  files  from  the 
Natural  Heritage  Programs  of  six  of  the  seven  states  where  L.  subcoriacea  is 
known  to  occur.  Unfortunately,  Louisiana  was  unable  to  release  the  location 
of  the  single  population  within  its  boundaries  due  to  a  limited  disclosure 
agreement  with  the  property  owner.  Collectively,  126  extant  L.  subcoriacea 
populations  were  available  for  potential  use  in  model  development.  These 
occurrence  data  generally  have  a  high  level  of  quality  control,  and  in  many 
cases  there  are  herbarium  records  associated  with  the  populations. 

2.1.2  Environmental  data 

For  each  900  m2  grid  cell  in  the  study  extent,  values  for  52  initial  environ¬ 
mental  predictors  representing  climate,  soils,  vegetation,  topography,  and 
fire  regime  were  calculated.  Climate  data  were  obtained  from  WorldClim 
version  1.4  (Hijmans  et  al.  2005).  This  dataset  includes  19  bioclimatic  in¬ 
dicators  that  represent  a  range  of  temperature  and  precipitation  summar¬ 
ies  (e.g.  trends,  seasonality,  and  extremes)  at  a  1  km  resolution  (appendix 
A).  These  data  were  resampled  using  a  simple  cell  assignment  to  match  the 
30  m  resolution  of  the  other  environmental  variables.  Sixteen  initial  soil 
variables  were  included  from  the  gSSURGO  (20 14)  soil  database  (Appen¬ 
dix  A),  a  rasterized  version  of  SSURGO  (USDA  20 15) .  Additional  vegeta¬ 
tion,  topographic,  and  fire  regime  variables  were  obtained  from  LAND- 
FIRE*  version  1.3  (USDOI  2013). 

2.1.3  MaxEnt  model  building 

The  maximum  entropy  method  was  used  as  implemented  in  MaxEnt  (Phil¬ 
lips,  Anderson,  and  Schapire  2006)  to  construct  predictive  models  of  L.  sub¬ 
coriacea  relative  probabilities  of  presence  (interpreted  as  habitat  suitability) 
across  the  southeastern  United  States.  Maximum  entropy  is  a  machine 
learning  approach  to  constructing  predictive  models  that  attempts  to  max¬ 
imize  the  entropy,  or  uncertainty,  of  the  probability  distribution  estimated 
from  the  data.  The  maximum  entropy  method  has  numerous  advantages, 
not  the  least  of  which  is  that  it  allows  the  modeling  of  presence- only  data.  At 


*  Landscape  Fire  and  Resource  Management  Planning  Tools  (LANDFIRE) 
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the  present  time,  maximum  entropy  is  the  most  popular  method  for  SDM 
and  performs  well  relative  to  other  techniques  (Elith  et  al.  2006) . 

Because  no  independent  testing  datasets  were  available,  5- fold  cross  valida¬ 
tion  was  used  to  partition  the  dataset  into  training  and  testing  groups  using 
ENMeval  version  0.21  (Muscarella  et  al.  2014)  within  R  (R  CoreTeam 
20 13) .  The  area  under  the  curve  (AUC)  of  a  receiver  operating  curve  (ROC) 
was  used  as  a  threshold- independent  metric  for  assessing  model  perfor¬ 
mance.  AUC  values  range  from  0- 1,  and  represent  the  probability  that  the 
model  ranks  a  random  presence  location  higher  than  a  random  background 
site  (Phillips,  Anderson,  and  Schapire  2006;  Phillips  and  Dudlk  2008). 
Overall,  the  AUC  metric  provides  an  assessment  of  how  accurately  the 
model  predicts  the  probability  of  occurrence  for  a  species  within  a  given 
area.  Models  with  AUC  values  greater  than  0.75  may  be  useful  for  identify¬ 
ing  the  potential  distribution  of  a  species  (Swets  1988,  Elith  et  al.  2006). 

Potential  sampling  bias  (Araujo  and  Guisan  2006,  Wintle  and  Bardos  2006, 
Veloz  2009,  Hijmans  2012,  Syfert,  Smith,  and  Coomes  2013),  number  of 
pseudoabsences  (McPherson,  J  etz,  and  Rogers,  2004;  Mateo  et  al.  2010, 
Iturbide  et  al.  2015),  and  background  extent  (VanDerWal  et  al.  2009,  Bar- 
bet- Massin  et  al.  2012)  are  all  known  to  influence  model  performance.  Thus, 
the  initial  modeling  of  the  training  data  and  full  set  of  52  environmental 
predictors  consisted  of  assessing  the  effects  of  altering  the  minimum  dis¬ 
tance  between  presence  observations,  pseudoabsences  and  the  sampling  ex¬ 
tent  on  the  AUC  values.  Spatially  thinned  datasets  were  created  with  pres¬ 
ence  and  randomly- placed  pseudoabsence  points  thinned  to  2,  5,  7,  and  10 
km  from  the  nearest  neighbor  using  the  R  package  spThin  (Aiello- Lammens 
et  al.  20 15) .  The  number  of  pseudoabsences  were  limited  to  twice  as  many 
as  the  presence  points.  The  effect  of  background  extent  on  AUC  at  different 
distances  ranging  from  2  to  200  km  from  presences  was  assessed.  The  back¬ 
ground  extent  of  pseudoabsence  point  selection  was  determined  by  observ¬ 
ing  the  distance  at  which  increases  in  AUC  were  minimal  ( VanDerW al  et  al. 
2009,  Iturbide  et  al.  20 15) .  This  process  identified  an  optimal  thinning  dis¬ 
tance  of  10  km  and  a  background  extent  encompassing  the  area  represent¬ 
ing  100  km  from  presences  (Figure  1) .  After  identifying  an  appropriate  thin¬ 
ning  distance,  sampling  extent,  and  number  of  pseudoabsences  with  the  full 
set  of  environmental  predictors,  the  full  set  of  environmental  predictors  was 
reduced  to  avoid  overfitting.  A  stepwise  process  was  used  in  which  predic¬ 
tors  that  had  a  <1%  mean  contribution  and  Pearson's  correlation  coeffi¬ 
cients  >|  0.9|  were  excluded,  and  the  correlated  predictor  with  the  highest 
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contribution  score  was  retained  (J  ueterbock  et  al.  20 15) .  The  process  was 
continued  until  the  contribution  scores  of  all  predictors  were  greater  than 
1%  and  the  mean  AUC  did  not  drop  below  the  arbitrary  AUC  threshold  of 
0.80.  In  addition,  model  overfitting  was  assessed  by  examining  the  differ¬ 
ence  between  training  and  test  data  AUC.  Warren  and  Seifert  (2011)  suggest 
that  overfit  models  should  generally  perform  well  on  training  data,  but 
poorly  on  test  data.  Conseguently,  when  the  training  and  test  data  AUC  ex¬ 
hibit  little  difference,  it  is  unlikely  one's  model  is  over- parameterized  such 
that  it  is  overly  specific  to  the  training  data. 


Figure  1.  Evaluation  of  the  influence  of  potential  sampling  bias  and  background  extent  on 

model  performance. 
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2.2  Model  application 

2.2.1  Areas  warranting  additional  survey  effort 

Several  approaches  for  identifying  locations  to  target  for  new  population 
discovery  have  been  proposed,  but  they  invariably  rely  on  identifying  some 
minimum  threshold  of  habitat  suitability.  One  approach  is  to  use  an  inde¬ 
pendent  dataset  to  identify  thresholds  of  suitability  associated  with  known 
occurrences  (e.g.,  Williams  et  al.  2009,  Crall  etal.  2013,  Wang  etal.  2014). 
However,  truly  independent  datasets  of  species  occurrences  are  rarely 
available;  even  occurrences  thinned  during  model  development  are  not 
suitable  for  this  intended  use.  Alternately,  threshold  values  have  been 
identified  by  examining  suitability  values  of  known  occurrences  used  for 
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model  development,  or  based  on  some  arbitrarily  proportion  (e.g.  top  dec¬ 
ile  or  quartile)  of  suitability  values  (e.g.,  Williams  et  al.  2009).  Because  in¬ 
dependent  data  on  L.  subcoriacea  occurrences  were  not  available,  the  up¬ 
per  quartile  (top  25%,  >0.73)  of  suitability  values  associated  with  occur¬ 
rences  used  for  model  development  to  identify  a  minimum  threshold.  Alt¬ 
hough  somewhat  arbitrary,  this  threshold  value  provided  a  balance  be¬ 
tween  over-  and  under- estimating  potential  survey  effort,  by  limiting  the 
area  and  locations  that  might  be  targeted  for  surveys. 

Additionally,  this  work  sought  to  place  emphasis  on  highly  suitable  habitat 
that  occurs  in  areas  that  have  not  been  subjected  to  extensive  anthropo¬ 
genic  degradation  and  fragmentation.  Therefore,  an  index  of  local  connect¬ 
edness  was  also  used  as  an  additional  criterion  to  identify  sites  warranting 
survey  effort.  Local  connectedness  characterizes  whether  the  connections 
within  natural  ecosystems  are  diminished  due  to  anthropogenic  barriers 
(i.e.  land  cover  conversion)  of  variable  permeability.  Specifically,  the  met¬ 
ric  used  was  one  developed  by  The  Nature  Conservancy  (TNC),  which  as¬ 
signs  resistance  weights  to  different  land  cover  types  based  on  the  degree 
to  which  natural  processes  have  likely  been  modified  (Anderson  et  al. 

20 16) .  The  weighting  scheme  assigns  low  resistance  to  any  adj  acent  natu¬ 
ral  cover  types  (e.g.,  forest,  wetland,  grassland),  and  higher  values  to  hu¬ 
man  altered  cover  types  (Table  1) .  Land  cover  data  used  to  develop  the 
metric  were  primarily  roads  data  from  the  20 14  Tiger  Road  dataset  (U.S. 
Census  2014)  and  the  2011  U.S.  Geological  Survey  (USGS)  National  Land 
Cover  Data  (NLCD),  which  represents  classified  Landsat  TM*  satellite  data 
for  the  48  conterminous  United  States  within  16  classes  (Homer  et  al. 

20 15) .  For  each  cell  of  the  resistance  grid,  the  values  of  neighboring  cells 
were  evaluated  out  to  a  distance  of  3  km  using  resistance  kernel  analysis 
(Compton  et  al.  2007) .  Additional  details  can  be  found  within  Anderson  et 
al.  ( 20 16) .  The  raster  layer  of  regionally  standardized  local  connectedness 
generated  by  TNC  using  this  analysis  had  a  grid  cell  size  of  90  m,  which 
was  resampled  to  match  the  30  m  resolution  of  the  habitat  suitability  map. 
ArcGIS  version  10.3.1  (ESRI,  Redlands,  CA,  USA)  was  used  to  identify  lo¬ 
cations  having  both  high  suitability  (>0.73)  and  connectedness  (>  average) 
as  potential  targets  for  new  population  discovery. 


*  Thematic  Mapper 
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Table  1.  Land  cover  classes  and  resistance  weights  used  to  estimate  local  connectedness. 


Land  cover  description 

Resistance  score 

Source 

Developed,  Open  Space 

8 

NLCD  2011 

Developed,  Low  intensity 

8 

NLCD  2011 

Developed,  Medium  Intensity 

9 

NLCD  2011 

Developed,  High  Intensity 

20 

NLCD  2011 

Barren  Land,  non-natural 

9 

NLCD  2011 

Herbaceous 

1 

NLCD  2011 

Hay/Pasture  (Coastal  Plain  &  Piedmont) 

3 

NLCD  2011 

Hay/Pasture  (Mountains) 

5 

NLCD  2011 

Cultivated  Crops 

7 

NLCD  2011 

Open  Water,  Shoreline  Distance  200-400m 

3 

NLCD  2011 

Open  Water,  Shoreline  Distance  >400  meters 

5 

NLCD  2011 

Major  Roads 

20 

Tiger  2014  (U.S.)  &  Open 
Street  Map  2014  (CA) 

Minor  Roads 

10 

Tiger  2014  (U.S.)  &  Open 
Street  Map  2014  (CA) 

Dirt  Roads 

Resistance  +1 

Open  Street  Map  2014 

Transmission  Lines 

9 

Ventex  2014 

Pipelines 

9 

Ventex  2014 

Railroads 

9 

CTS  2015 

Unprotected/Private  Industrial  Forest  (U.S.) 

3 

SE  GAP,  Parcelpoint,  OSI* 

Protected  Industrial  Forest  (U.S.) 

1.5 

SE  GAP,  Parcelpoint,  OSI 

‘Open  Systems  Interconnection  (OSI) 

NLCD  2011 


NLCD  2011 


NLCD  2011 


Woody  Wetlands 


Emergent  Herbaceous  Wetlands 


Open  Water,  Shoreline  Distance  <200  m 


NLCD  2011 

NLCD  2011 

NLCD  2011 

NLCD  2011 

NLCD  2011 

Barren  Land,  natural 

Deciduous  Forest 

Evergreen  Forest 

Mixed  Forest 

Shrub/Scrub 
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2.2.2  Areas  for  potential  conservation  and  management 

Three  types  of  conservation  and  management  efforts  were  identified  as 
likely  to  be  of  interest  to  public  and  private  entities  concerned  about  con¬ 
serving  known  L.  subcoriacea  occurrences.  The  first  two  are  related  to 
identifying  sites  for  either  acguisition  or  development  of  conservation 
agreements,  while  the  third  addresses  identification  of  occupied  sites  po¬ 
tentially  needing  habitat  restoration  and/  or  population  monitoring.  Occu¬ 
pied  sites  potentially  suitable  for  acguisition  or  establishing  conservation 
agreements  include  those  that  are  of  high  suitability,  unprotected,  and 
having  average  to  high  connectedness  of  natural  land  cover  types.  Unoccu¬ 
pied,  unprotected  sites,  having  high  suitability  and  average  to  high  con¬ 
nectedness,  but  near  known  occurrences  (<1  km),  can  potentially  act  as 
conservation  buffers  and  may  also  be  suitable  targets  for  acguisition  or  es¬ 
tablishing  conservation  agreements.  Occupied  sites  potentially  needing 
habitat  restoration  and/  or  population  monitoring,  are  those  that  are  pro¬ 
tected  as  a  first  priority,  having  average  to  high  connectedness,  but  in  loca¬ 
tions  modeled  to  have  moderate  to  low  suitability.  ArcGIS  was  used  to 
identify  locations  meeting  these  criteria.  Information  on  land  protections 
was  sourced  from  the  USGS,  Gap  Analysis  Program,  Protected  Areas  Data¬ 
base,  version  1.4  (http://gapanaivsis.usgs.gov/padus/) .  Data  for  the  other  criteria 
have  been  described  previously. 

2.2.3  Areas  suitable  for  potential  (re)introduction 

In  addition  to  managing  extant  populations,  (re)introduction  can  form  an 
important  part  of  an  integrated  strategy  for  rare  plant  conservation  (e.g., 
Menges  and  Kennedy  2007,  Thorpe  and  Kaye  2011,  Halsey  et  al  2015). 
However,  lack  of  knowledge  about  the  suitability  of  potential  sites  for 
(re)introduction  can  either  hinder  the  willingness  of  land  managers  and 
regulators  to  adopt  the  approach,  or  lower  success  rates  (Godefroid  et  al. 
2011,  Drayton  and  Primack  2012,  Questad  et  al.  2014).  Therefore,  this 
work  sought  to  identify  sites  potentially  suitable  for  (re)introduction  to 
promote  the  viability  of  this  conservation  strategy.  Sites  potentially  suita¬ 
ble  for  (re)introduction  include  those  that  are  unoccupied,  highly  suitable, 
protected,  and  having  average  to  high  connectedness.  ArcGIS  was  used  to 
identify  locations  meeting  these  criteria. 
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3  Results  and  Discussion 

3.1  Model  performance 

The  final  model  provided  "good"  differentiation  between  occupied  and 
background  locations  (Swets  1988);  the  AUC  egualed  0.80.  The  test  of 
model  overfitting  also  yielded  satisfactory  results,  as  the  training  and  test¬ 
ing  AUC  were  comparable  (AUC  train  =0.84  ±  0.01,  AUC  test  =  0.80  ±0.10). 
Twenty  environmental  predictors  representing  all  five  predictor  types  were 
retained  in  the  final  model  after  predictor  reduction  (Appendix  A) .  Among 
the  five  types  of  environmental  predictors,  vegetation  predictors  only  ac¬ 
counted  for  four  of  the  final  20,  but  collectively  contributed  42%,  while  fire 
regime  predictors  were  the  second- most  important  category,  collectively 
contributing  23%  (Table  2).  The  other  three  types  of  environmental  pre¬ 
dictors  (climate,  soil,  and  topography)  had  lower  percent  contributions 
that  ranged  from  5- 15%. 

Table  2.  Five  predictor  types,  number  of  predictors  retained,  and 
their  relative  importance  in  the  final  L.  subcoriacea  model. 


Predictor  type 

Number  retained 

Percent  contribution 

Vegetation 

4 

42.3 

Fire  regime 

4 

23.3 

Soils 

4 

14.9 

Climate 

5 

14.0 

Topography 

3 

5.6 

Considering  individual  environmental  predictors,  biophysical  settings 
(BPS)  and  vegetation  condition  class  (VCC)  emerged  as  the  most  im¬ 
portant.  All  other  predictors  contributed  <10%  in  the  final  model.  BPS  is  a 
vegetation  related  predictor  while  VCC  is  classified  as  a  fire  predictor  (Ap¬ 
pendix  A) .  BPS  was  by  far  the  most  influential  (31.5%  contribution),  con¬ 
tributing  over  twice  as  much  as  the  second  most  important  environmental 
predictor,  VCC.  BPS  represents  the  vegetation  that  may  have  been  domi¬ 
nant  on  the  landscape  before  European  settlement  and  is  based  on  both 
the  current  biophysical  environment  and  an  approximation  of  the  histori¬ 
cal  disturbance  regime  (USDOI  2013).  Locations  occupied  by  L.  subcoria¬ 
cea  were  represented  by  14  different  BPS  categories,  but  only  four  catego¬ 
ries  accounted  for  >10%  of  known  occurrences  (Table  3) . 
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Table  3.  Biophysical  settings  of  L  subcoriacea  populations  and  their  proportional 

representation. 


Biophysical  settings 

Percent  of  L. 
subcoriacea 

occurrences 

Atlantic  Coastal  Plain  Fall-line  Sandhills  Longleaf  Pine  Woodland 

30.3 

East  Gulf  Coastal  Plain  Interior  Upland  Longleaf  Pine  Woodland 

15.9 

Gulf  and  Atlantic  Coastal  Plain  Small  Stream  Riparian  Systems 

13.6 

Atlantic  Coastal  Plain  Streamhead  Seepage  Swamp-Pocosin-Baygall 

12.1 

Atlantic  Coastal  Plain  Dry  and  Dry-Mesic  Oak  Forest 

6.8 

Atlantic  Coastal  Plain  Mesic  Hardwood  Forest 

6.8 

Southern  Coastal  Plain  Seepage  Swamp  and  Baygall 

3.0 

Atlantic  Coastal  Plain  Upland  Longleaf  Pine  Woodland 

2.3 

East  Gulf  Coastal  Plain  Savanna  and  Wet  Prairie 

2.3 

Southern  Piedmont  Dry  Oak(-Pine)  Forest 

2.3 

Central  Atlantic  Coastal  Plain  Wet  Longleaf  Pine  Savanna  and  Flatwoods 

1.5 

Gulf  and  Atlantic  Coastal  Plain  Swamp  Systems 

1.5 

Florida  Longleaf  Pine  Sandhill 

0.8 

Southern  Coastal  Plain  Mesic  Slope  Forest 

0.8 

Additional  vegetation- related  predictors  retained  in  the  final  model  in¬ 
cluded  environmental  site  potential  (ESP),  canopy  height  (CH),  and  canopy 
cover  (CC).  ESP  represents  the  vegetation  that  could  become  established 
during  late  or  climax  stages  of  successional  development  at  a  given  site 
based  on  the  biophysical  environment  (USDOI  2013).  Locations  occupied 
by  L.  subcoriacea  represented  12  ESP  categories  (Table  4),  which  largely 
paralleled  BPS  categories;  only  the  same  four  categories  accounted  for  >10% 
of  known  occurrences.  CC  and  CH  at  locations  occupied  by  L.  subcoriacea 
were  highly  variable,  ranging  from  0-100%  (median  =60-70%)  and  0-50  m 
(median  =  10-25  m),  respectively.  CC  and  CH  influence  light  availability 
(Martens,  Breshears,  and  Meyer  2000),  which  can  affect  plant  photosyn¬ 
thetic  rates,  growth,  herbivory,  water  balance,  and  reproduction.  The  effects 
of  light  availability  are  well- documented  for  L.  benzoin  (Niesenbaum  1992, 
1993;  Luken  et  al.  1997;  Muth  et  al.  2008),  and  light  availability  has  also 
been  suggested  to  affect  L.  subcoriacea  plant  vigor  and  reproduction  (pers. 
obs.,  Tom  Patrick) .  CC  and  CH  not  only  influence  light  availability,  but  also 
fuels  and  microclimate,  which  have  been  shown  to  affect  fire  return  interval, 
behavior  and  fire- vegetation  feedbacks  in  communities  occupied  by  L.  sub¬ 
coriacea  (J  ust,  Hohmann,  and  Hoffmann  2016). 
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Vegetation- related  predictors  at  locations  occupied  by  L.  subcoriacea  match 
field- based  descriptions  of  the  species'  known  associations  with  perma¬ 
nently  moist  to  wet,  evergreen  shrub- dominated  seepage  wetlands,  namely 
Atlantic  white  cedar  forests,  coastal  plain  small  stream  swamps,  streamhead 
pocosins,  and  sandhill  seeps  in  North  and  South  Carolina  (Some,  Gray,  and 
Crutchfield  2006,  Wall  et  al.  2013),  and  the  wettest  portions  of  sphagnous 
bogs  in  Mississippi  and  Florida  (Gordon,  J  ones,  and  Wiseman  1986) . 


Table  4.  Environment  site  potential  categories  of  L.  sucoriacea  populations  and  their 

proportional  representation. 


Environment  site  potential 

Percent  of  L. 
subcoriacea 

occurrences 

Atlantic  Coastal  Plain  Fall-line  Sandhills  Longleaf  Pine  Woodland 

25.6 

East  Gulf  Coastal  Plain  Interior  Upland  Longleaf  Pine  Woodland 

16.8 

Atlantic  Coastal  Plain  Streamhead  Seepage  Swamp-Pocosin-Baygall 

16.0 

Gulf  and  Atlantic  Coastal  Plain  Small  Stream  Riparian  Systems 

15.2 

Atlantic  Coastal  Plain  Dry  and  Dry-Mesic  Oak  Forest 

8.8 

Atlantic  Coastal  Plain  Mesic  Hardwood  Forest 

7.2 

Southern  Coastal  Plain  Seepage  Swamp  and  Baygall 

3.2 

Gulf  and  Atlantic  Coastal  Plain  Swamp  Systems 

2.4 

Southern  Coastal  Plain  Mesic  Slope  Forest 

1.6 

Southern  Piedmont  Mesic  Forest 

1.6 

Central  Atlantic  Coastal  Plain  Wet  Longleaf  Pine  Savanna  and  Flatwoods 

0.8 

Florida  Longleaf  Pine  Sandhill 

0.8 

The  most  important  fire  regime  predictor  was  VCC  (12.5%  contribution) . 
VCC,  which  represents  a  simple  categorization  of  vegetation  departure 
(VDEP),  and  denotes  the  degree  to  which  current  vegetation  differs  from 
historical  vegetation  reference  conditions  (USDOI  2013).  VDEP  at  loca¬ 
tions  occupied  by  L.  subcoriacea  ranged  from  38-81%  (mean  =  50.5%,  sd 
=  ±7.9)  Additional  fire  regime  predictors  included  fire  regime  groups 
(FRG)  and  percentage  of  replacement  severity  fires  (PRS)  (Appendix  A). 
FRG  characterizes  the  presumed  historical  fire  regimes  within  landscapes 
based  on  interactions  between  vegetation  dynamics,  fire  spread,  fire  ef¬ 
fects,  and  spatial  context  (USDOI  2013).  Locations  occupied  by  L.  subcori¬ 
acea  include  three  different  FRG  categories  (FRG  I,  III  and  V),  which  span 
fire  return  intervals  of  <35  to  >200  years.  However,  the  most  prevalent 
category  was  FRG  I,  which  represents  a  <  35-year  fire  return  interval,  with 
low  and  mixed  severity  fire.  PRS  ranged  from  0- 15%  at  locations  occupied 
by  L.  subcoriacea  (median  =  0-5%),  with  replacement  severity  defined  as 
greater  than  75%  average  top-kill  within  a  typical  fire  perimeter  for  a  given 
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vegetation  type.  Although  caution  should  be  taken  not  to  infer  causality 
from  the  correlative  approach  used  by  MaxEnt  to  predict  habitat  suitabil¬ 
ity,  it  is  worth  noting  that  the  fire  variables  retained  in  the  final  model  and 
their  values  match  those  reported  for  L.  subcoriacea  (Some,  Gray,  and 
Crutchfield  2006,  Wall  et  al.  2013).  Fire  history  is  a  considerable  determi¬ 
nant  of  the  dominant  vegetation  of  wetlands  occupied  by  L.  subcoriacea 
(Some,  Gray,  and  Crutchfield  2006),  which  are  estimated  to  have  a  fire  re¬ 
turn  interval  of  7-50  years  (Frost  1998,  Stambaugh,  Guyette,  and  Mar- 
schall  2011,  J  ust,  Hohmann,  and  Hoffmann  2016). 

Soil  predictors  retained  in  the  final  model  included  field  capacity  (WTENTH- 
BAR),  depth  to  wet  soil  layer  (WTDEPANNMIN),  soil  permeability  (KSAT), 
and  hydrogen  ion  activity  of  soil  (PH1T01H20) .  The  first  three  of  these  pre¬ 
dictors  are  related  to  soil  moisture,  while  the  last  reflects  soil  pH  (Table  5) . 
Soil  pH  affects  many  aspects  of  plant  physiology,  including  nutrient  balance 
and  growth  rates.  PH1T01H20  made  the  largest  contribution  at  7.58%,  fol¬ 
lowed  by  KSAT  (3.68%),  WTDEPANNMIN  (2.17%),  and  WTENTHBAR 
( 1.42%) .  Locations  occupied  by  L.  subcoriacea  generally  had  higher  hydraulic 
conductivity,  and  a  lower  depth  to  the  water  table  (WTDEPANMIN)  relative 
to  unoccupied  locations,  though  the  differences  were  not  significant.  These 
values  match  field  descriptions  of  acidic,  permanently  saturated,  high  organic 
matter  soils  (Gordon,  J  ones,  and  Wiseman  1986). 

Climate- related  predictors  retained  in  the  final  model  included  isother- 
mality  (BI03),  mean  temperature  of  driest  guarter  (BI09),  annual  precipi¬ 
tation  (BI012),  precipitation  seasonality  (BI015),  and  precipitation  of  the 
driest  guarter  (BI017) .  As  might  be  anticipated  by  the  species'  geographic 
range,  the  dimate  of  occupied  locations  was  moist  temperate,  with  modest 
seasonal  variation  in  temperatures  and  predpitation  (Table  5) . 


Table  5.  Mean,  standard  deviation  and  range  of  continuous  soil,  topographic  and  climate 
environmental  predictors  at  locations  where  L.  subcoriacea  populations  occur. 


Variable 

Units 

Mean 

Standard 

deviation 

Range 

ASPECT 

None 

0.50 

0.72 

-1-1 

ESI012 

mm 

1303 

184.6 

1142-1645 

BI015 

None 

18.9 

2.2 

12-23 

ESI017 

mm 

258 

38.8 

215-331 

BI03 

% 

40 

1.74 

38-44 

BI09 

°c 

13.8 

3.4 

10.4-22.9 

KSAT 

pm/s 

889.6 

392.2 

0-1283 
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Variable 

Units 

Mean 

Standard 

deviation 

Range 

PH1T01H20 

None 

28.8 

13.4 

0-88.5 

SLOPE 

% 

2.63 

2.49 

0-13 

TPI 

None 

-0.44 

1.24 

-3.48  -  4.04 

WTDEPANNMIN 

cm 

11.7 

28.7 

0-122 

WTENTHBAR 

% 

27.3 

15.1 

0-46.8 

Topographic  predictors  retained  in  the  final  model  included  slope,  aspect, 
and  topographic  position  index  (TPI).  Locations  occupied  by  L.  subcoria- 
cea  were  predominantly  mid-positon,  shallow  slopes  (Table  5),  which 
matches  known  associations  with  mid- slope,  hillside  seepage  wetlands  in 
the  sandhills  region  of  North  and  South  Carolina  (Some,  Gray,  and 
Crutchfield  2006,  Wall  et  al.  2013),  and  level  to  slightly  sloping  terrain 
within  the  Gulf  Coastal  Plain  of  Alabama,  Mississippi,  and  Louisiana.  Alt¬ 
hough  considerable  variation  in  the  aspect  of  occupied  sites  was  found,  the 
majority  were  north- facing  (Table  5) .  Approximately  60%  of  sites  had  val¬ 
ues  close  to  1,  which  corresponds  to  north,  while  only  20%  were  either 
south- facing,  or  east- west  facing.  This  association  with  north- facing  slopes 
is  a  novel  insight  about  site  suitability. 

3.2  Habitat  suitability  maps 

Potential  habitat  with  high  suitability  (>0.73)  was  distributed  throughout 
the  region  (Appendix  B),  with  2%  of  the  region  identified  as  highly  suitable. 
Percent  area  of  high  suitability  ranged  from  10%  in  the  Northeastern  Pied¬ 
mont  physiographic  region  to  less  than  .0 1%  in  the  Mississippi  Delta  and 
the  Southeastern  Coastal  Plain  physiographic  regions  (Appendix  B).  The 
high  values  for  the  Northeastern  Piedmont  physiographic  region  are  inter¬ 
esting;  however  most  of  the  area  is  outside  the  known  historical  range  of  L. 
subcoriacea,  with  only  a  few  known  locations  at  the  southeastern  edge  of 
the  ecoregion.  North  Carolina  had  the  highest  percentage  of  suitable  habitat 
(5.0%)  and  Georgia  had  the  lowest  ( <0.01%).  Mississippi  had  the  highest 
percentage  (24.0%)  of  high  suitability  area  on  protected  lands,  while  Louisi¬ 
ana  and  North  Carolina  had  the  lowest  (4.6  and  4.9%,  respectively). 

Considering  DoD  installations  known  to  support  L.  subcoriacea  popula¬ 
tions  and  the  landscapes  within  a  10  km  buffer  of  the  installation  bounda¬ 
ries,  Eglin  AFB  and  Camp  Shelby  J  FTC  were  found  to  have  large  areas  of 
highly  suitable  habitat  compared  to  Fort  Bragg  and  Fort  J  ackson  (Tables  6 
and  7) .  Except  for  Camp  Shelby  J  FTC,  much  of  the  highly  suitable  habitat 
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within  the  buffers  of  installations  is  unprotected.  However,  there  are  small 
amounts  of  state- protected  lands  with  high  suitability  near  Fort  Bragg  and 
Eglin  AFB  (Table  7) . 

Figure  2.  Ecoregions  occupied  by  Lindera  subcoriacea  in  the  southeastern  United  States. 


Table  6.  Area  (Ha)  on  installations  representing  highly  suitable  L.  subcoriacea  habitat, 
locations  warranting  additional  survey  effort,  and  sites  potentially  suitable  for  (re)introduction. 


Installation 

Total  area 

High  suitability 

Additional  survey 

(Re)introduction 

Eglin  AFB 

188041.1 

2644.5 

1900.3 

1888.9 

Fort  Bragg 

65941.0 

139.1 

41.0 

37.6 

Camp  Shelby 

55995.9 

1859.0 

1595.1 

1489.7 

Fort  Jackson 

21004.9 

20.3 

10.4 

10.0 
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Table  7.  Areas  (Ha)  within  a  10  km  buffer  of  installation  boundaries  representing  highly 
suitable  L.  subcoriacea  habitat,  locations  warranting  additional  survey  effort  and  sites 
potentially  suitable  for  (re)introduction.  Values  are  summarized  by  type  of  protection. 


Installation 

Protection 

Type 

Total  Area 

High 

Suitability 

Additional 

Survey 

(Re) 

Introduction 

Camp  Shelby 

Federal 

36871.11 

3237.75 

3187.8 

3175.11 

NGO 

1952.19 

210.78 

210.78 

210.51 

None 

108747.45 

5081.76 

4110.93 

0 

Private 

1018.80 

1.8 

1.8 

1.8 

State 

6073.56 

379.08 

293.04 

292.59 

Unknown 

4.95 

0.18 

0.18 

0 

total  = 

154668.06 

8911.35 

7804.53 

3680.01 

Eglin  AFB 

Federal 

4710.33 

1.35 

0 

0 

Local 

144.90 

0.09 

0 

0 

NGO 

13.41 

0 

0 

0 

None 

187903.62 

4535.01 

1095.12 

0 

Private 

7418.34 

1.53 

0.63 

0.63 

State 

34575.03 

282.24 

98.19 

97.11 

Unknown 

661.14 

0.72 

0.45 

0.45 

total  = 

235426.77 

4820.94 

1194.39 

98.19 

Fort  Bragg 

Federal 

156.69 

5.58 

5.58 

5.58 

Local 

580.95 

1.26 

0 

0 

NGO 

630.54 

1.53 

0 

0 

None 

164145.51 

681.48 

262.08 

0 

Private 

1377.54 

8.28 

3.15 

2.97 

State 

19083.96 

380.07 

293.58 

258.39 

Unknown 

334.89 

0.9 

0.45 

0.36 

total  = 

186310.08 

1079.1 

564.84 

267.3 

Fort  Jackson 

Local 

507.42 

3.15 

0 

0 

None 

86765.94 

212.13 

37.98 

0 

Private 

2296.08 

5.31 

5.31 

5.31 

State 

949.05 

4.68 

0 

0 

Unknown 

798.21 

0 

0 

0 

total  = 

91316.70 

225.27 

43.29 

5.31 
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3.3  Model  application 

3.3.1  Areas  warranting  additional  survey  effort 

A  total  of  4719.3  km2  was  identified  as  warranting  additional  survey  effort. 
These  areas  were  variably  distributed  across  seven  states,  seven  ecore- 
gions,  and  five  consolidated  types  of  land  managers  (Appendix  B) .  North 
Carolina  had  the  most  area  ( 1522  km2,  33%  of  the  total)  identified  for  ad¬ 
ditional  survey  effort,  which  is  almost  entirely  located  on  unprotected  pri¬ 
vate  lands  within  the  Northeastern  Piedmont  ecoregion  ( 1482  km2) .  Ala¬ 
bama  and  Mississippi  also  had  sizable  areas  identified  for  additional  sur¬ 
vey  effort,  1119  and  1077  km2,  respectively.  Georgia  had  the  least  amount 
of  area  (99.5  km2)  identified  for  additional  survey  effort  and  77%  of  this 
area  is  on  private  lands.  Among  ecoregions,  the  largest  areas  and  93%  of 
the  total  identified  area  warranting  additional  survey  effort  were  located  in 
the  Coastal  Gulf  Plains  (1603  km2).  Northeastern  Piedmont  (1567  km2), 
and  Gulf  Plains  ( 1237  km2) . 

Considering  DoD  installations  known  to  support  L.  subcoriacea  popula¬ 
tions  and  the  landscapes  within  a  10  km  buffer  of  the  installation  bounda¬ 
ries,  Camp  Shelby  J  FTC  and  Eglin  AFB  were  found  to  have  relatively  large 
areas  warranting  additional  survey  effort  both  on  the  installation  and 
within  the  surrounding  landscape  (Tables  6  and  7).  In  contrast.  Fort  Bragg 
and  Fort  J  ackson  had  much  smaller  areas  within  the  installation  bounda¬ 
ries  and  the  surrounding  landscapes.  For  example,  one  can  see  widely  dis¬ 
tributed  areas  on  Camp  Shelby  J  FTC  and  within  the  neighboring  DeSoto 
National  Forest  where  additional  searches  for  L.  subcoriacea  may  be  suc¬ 
cessful  (Figures  3  and  4). 
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Figure  3.  Example  locations  on  and  near  Camp  Shelby  JFTC  potently  warranting  additional 

survey  effort. 
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Figure  4.  Close  up  of  locations  on  Desoto  National  Forest  near  Camp  Shelby  JFTC  warranting 
additional  survey  effort  for  L.  subcoriacea. 
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3.3.2  Areas  for  potential  conservation  and  management 

These  analyses  focused  on  identifying  two  types  of  areas  of  potential  inter¬ 
est  for  targeted  acguisition  and/  or  development  of  conservation  agree¬ 
ments.  Unfortunately,  no  populations  under  private  ownership  within 
highly  suitable  and  connected  habitat  were  identified.  There  are  only  32 
known  L.  subcoriacea  populations  in  private,  unprotected  ownership;  all 
of  these  appear  to  be  in  degraded  habitat,  making  them  relatively  unat¬ 
tractive  for  acguisition  or  inclusion  within  a  conservation  agreement  with¬ 
out  consideration  of  their  surrounding  landscapes.  Consequently,  the 
landscapes  within  a  1  km  radius  buffer  of  these  unprotected  populations 
were  explored  to  identify  the  area  of  highly  suitable  and  connected  habitat 
that  might  facilitate  local  metapopulation  dynamics.  For  example,  lands 
near  unprotected  L.  subcoriacea  populations  in  Mississippi  that  could  be 
targeted  for  acquisition  or  development  of  conservation  agreements  can  be 
seen  in  (Figure  5).  The  areas  within  buffers  varied  across  individual  popu¬ 
lations;  Table  8  lists  the  number  of  populations  and  mean  area  by  state 
that  are  potentially  suitable  for  landscape- based  conservation  strategies. 
No  privately  owned,  connected,  high  suitability  habitat  near  known  L.  sub¬ 
coriacea  occurrences  were  identified  on  public  lands  as  possible  targets  for 
extending  existing  conservation  buffers. 

Table  8.  Numbers  of  L.  subcoriacea  populations  on  unprotected  lands  and  the  mean  area  of 
highly  suitable  and  connected  habitat  within  surrounding  1  km  radius  landscapes. 


State 

Number  of  populations 

Mean  area  (Ha) 

AL 

3 

4.2 

GA 

3 

5.0 

MS 

9 

5.1 

NC 

13 

6.0 

sc 

4 

2.7 
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Figure  5.  Example  of  locations  within  1  km  buffers  of  unprotected  L.  subcoriacea  populations 
in  Mississippi  that  may  be  suitable  for  acquisition  or  development  of  conservation 
agreements  as  part  of  landscape-based  conservation  strategies. 
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This  work  also  sought  to  identify  L.  subcoriacea  populations  on  protected 
lands  potentially  needing  habitat  restoration  or  population  monitoring.  Re¬ 
sults  of  these  analyses  suggest  that  46  of  88  (52%)  populations  on  protected 
lands  occur  in  locations  having  low  to  moderate  habitat  suitability,  but  high 
connectedness  (Table  9).  An  example  of  populations  on  Federal  lands  po¬ 
tentially  needing  habitat  restoration  or  population  monitoring  are  those  on 
the  U.S.  Department  of  Energy  (DOE)  Savannah  River  Site.  The  percent  of 
known  L.  subcoriacea  occurrences  identified  as  potentially  needing  restora¬ 
tion,  was  roughly  equal  among  Federal-  and  state- managed  land,  which  are 
responsible  for  the  vast  majority  of  populations.  Although,  one  can  alter¬ 
nately  interpret  these  findings  as  poor  model  performance,  it  is  generally 
known  that  the  condition  of  L.  subcoriacea  habitats  and  populations  can 
wax  and  wane  in  response  to  fire,  flooding,  and  other  environmental  varia¬ 
bles.  Consequently,  land  managers  should  consider  evaluating  the  condition 
of  populations  in  the  field  and  investing  effort  in  restoring  them  as  neces¬ 
sary  (e.g.,  via  prescribed  fire  or  thinning  the  canopy) . 


Table  9.  Number  and  percent  of  L.  subcoriacea  populations 
identified  as  potentially  needing  restoration  or  population 
monitoring. 


Land  manager 

Number  of  populations 

Percent  of  populations 

Federal 

73 

52 

State 

12 

50 

NGO 

2 

100 

Local 

1 

0 

These  results  suggest  that  modeling  habitat  suitability  can  provide  insight 
into  the  appropriateness  of  different  conservation  and  management  strate¬ 
gies  for  L.  subcoriacea  populations  on  protected  lands  versus  those  on  un¬ 
protected  lands. 

3.3.3  Areas  suitable  for  potential  (re)introduction 

A  total  of  906.9  km2  was  identified  as  potentially  suitable  for  (re)introduc- 
tion  efforts.  These  areas  were  variably  distributed  across  seven  states, 
seven  ecoregions,  and  five  consolidated  types  of  protected  lands  (Appendix 
B) .  Among  types  of  protected  lands,  the  largest  area  suitable  for  potential 
(re)introduction  is  located  on  Federal  properties  (477.7  km2).  Mississippi 
had  the  most  area  (344.7  km2)  identified  for  potential  (re)introduction, 
with  a  large  region  of  state- owned  land  (Pascagoula  State  Fish  and  Wildlife 
Area)  concentrated  in  the  Coastal  Gulf  Plains  ecoregion  (Figure  6).  Among 
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ecoregions,  the  Coastal  Gulf  Plains  had  the  largest  (371  km2)  and  the  Mis¬ 
sissippi  Delta  had  the  smallest  (0.0  km2)  areas  potentially  suitable  for 
(re)introduction. 

Considering  DoD  installations  known  to  support  L.  subcoriacea  popula¬ 
tions  (Table  6),  both  Camp  Shelby  J  FTC  and  Eglin  AFB  were  identified  to 
have  sizable  areas  suitable  for  (re)introduction.  In  contrast.  Fort  Bragg 
and  Fort  J  ackson  had  relatively  little  area  suitable  for  (re)  introduction. 
Considering  10  km  buffers  of  these  installations,  only  Camp  Shelby  J  FTC 
was  identified  to  have  a  sizeable  area  potentially  suitable  for  (re)introduc- 
tion  (Table  7),  and  these  areas  were  largely  under  Federal  ownership.  Alt¬ 
hough  the  available  areas  are  much  smaller,  there  were  also  some  sites  on 
state  game  lands  near  Fort  Bragg/  Camp  Mackall  that  could  contribute  to 
local  metapopulation  dynamics  if  (re)introductions  were  implemented 
(Figure  5). 

The  success  rates  of  rare  plant  (re)introductions  are  often  low  due  to  diffi¬ 
culty  identifying  suitable  habitat  for  restoration  (Godefroid  et  al.  20 11, 
Drayton  and  Primack  2012,  Questad  et  al.  2014).  Conseguently,  SDM  can 
play  an  important  role  in  guiding  listed  and  at-risk  plant  (re)introductions. 
Still,  there  are  relatively  few  examples  of  published  studies  where  MaxEnt 
has  been  used  to  identify  potential  (re)introduction  sites  for  listed  and  at- 
risk  plants  (e.g.,  Adhikari,  Barik,  and  Upadhaya  2012).  These  results  sug¬ 
gest  there  are  many  locations  potentially  suitable  for  (re)introducing  L. 
subcoriacea  across  its  range.  However,  in  states  (e.g.,  Georgia)  and  ecore¬ 
gions  where  suitable  (re)introduction  sites  are  less  available,  management 
options  may  be  limited  to  restoration  of  existing  occupied  sites. 

Two  hundred  L.  subcoriacea  seedlings  were  successfully  outplanted  into 
eight  (re)introduction  sites  on  Fort  Bragg.  Germination  rates  were  high 
(-85%)  and  1-year  (2016)  survival  rates  were  approximately  65%  across 
the  sites  (unpublished  data) .  The  relatively  high,  albeit  short-term,  survi¬ 
vorship  of  these  outplants  suggests  that  (re)introduction  could  prove  suc¬ 
cessful  for  either  augmenting  small  populations,  or  establishing  new  popu¬ 
lations  to  increase  the  metapopulation  structure  within  clusters  of  popula¬ 
tions.  These  outplanting  data  also  speak  to  the  potential  influence  of  seed, 
dispersal,  and  recruitment  limitation  on  the  realized  distribution  of  the 
species.  Experiments  were  initiated  to  robustly  test  the  differential  effects 
of  these  limitations  on  the  species'  population  dynamics. 
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Figure  6.  Large  area  potentially  suitable  for  L.  subcoriacea  (re)introduction  efforts  within  the 

Coastal  Gulf  Plains  ecoregion. 
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Figure  7.  Locations  on  state  game  lands  near  Fort  Bragg  and  Camp  Mackall  identified  as 
being  potentially  suitable  for  L.  subcoriacea  (re)introduction  efforts. 
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4  Conclusions  and  Recommendations 

4.1  Conclusions 

Species  distribution  modeling  can  be  a  valuable  tool  for  rare  plant  conserva¬ 
tion  assessment  and  planning.  This  application  of  SDM  to  L.  subcoriacea 
identified  abundant  high  suitability  habitat  across  the  species'  range,  numer¬ 
ous  locations  warranting  additional  survey  effort,  specific  conservation  as¬ 
sessment  and  management  needs  for  occupied  sites  and  their  surrounding 
landscapes,  and  sites  potentially  suitable  for  population  (re)introduction. 

Although  it  is  anticipated  that  these  results  will  benefit  multiple  facets  of 
L.  subcoriacea  conservation,  it  is  also  recognized  that  species  distributions 
are  affected  by  historical  processes,  dispersal  limitation,  and  the  dynamics 
of  natural  systems.  Collectively  these  factors  preclude  species  from  occu¬ 
pying  their  whole  potential  distribution  (Marcer  et  al.  20 13) .  This  can  gen¬ 
erate  mismatches  between  predicted  habitat  suitability  and  current  site  oc¬ 
cupation.  Therefore,  the  absence  of  L.  subcoriacea  from  areas  identified  as 
highly  suitable,  or  presences  in  areas  identified  as  unsuitable  cannot  be 
simply  interpreted  as  a  model  failure. 

4.2  Recommendations 

4.2.1  Potential  opportunities  for  model  improvement 

As  is  often  the  case,  several  opportunities  for  potential  model  improve¬ 
ment  came  to  light  after  completing  the  analyses.  First,  any  assessment  of 
habitat  suitability  should  be  considered  a  temporary  estimate  given  the  dy¬ 
namic  nature  of  the  vegetation  communities  occupied  by  L.  subcoriacea. 
The  species  occupies  habitats  that  are  in  constant  flux  due  to  the  influence 
of  fire,  fire  suppression,  and  flooding  (e.g.,  due  to  beaver  activity)  on  vege¬ 
tation  composition  and  structure.  This  effort  used  vegetation  and  disturb¬ 
ance  data  provided  within  the  most  recent  version  of  LANDFIRE  (version 
1.3),  which  represents  the  conditions  in  2012  and  the  landscape  disturb¬ 
ances  that  occurred  in  20 11  and  20 12.  An  updated  version  ( 1.4)  of  LAND- 
FIRE  is  scheduled  to  be  completed  in  March  2017  and  a  remapping  effort 
designed  to  produce  new  base  maps  of  the  LANDFIRE  product  suite  (ver¬ 
sion  2.0)  will  likely  be  completed  by  2020.  It  is  advised  that  the  model  be 
revised  when  these  newer  data  become  available  as  they  will  better  repre¬ 
sent  actual  conditions  in  the  field. 
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Second,  it  is  also  possible  that  some  of  the  locations  within  the  present  da¬ 
taset  used  for  model  development  are  no  longer  occupied.  Wall  et  al. 

(2013)  combined  multi-year  field  inventories  with  modeling,  and  esti¬ 
mated  that  only  72%  of  known  L.  subcoriacea  populations  have  persisted 
in  North  Carolina  over  the  past  30  years.  Fortunately  many  of  the  sites  in 
North  Carolina  and  all  known  sites  in  Georgia  were  recently  revisited 
(Wall  et  al.  2013;  pers.  com,  Tom  Patrick).  However,  last  observation  dates 
from  other  states  occasionally  went  as  far  back  as  1985  and  may  be  a 
source  of  modeling  error.  Therefore  a  separate  model  of  habitat  suitability 
generated  using  only  recent  presence  data  may  be  worth  pursuing. 

4.2.2  Implementing  SDM-informed  conservation 

Public  and  private  entities  concerned  about  L.  subcoriacea  conservation 
should  use  the  information  generated  by  this  effort  to  enhance  conserva¬ 
tion  assessment  and  planning  efforts.  The  USFWS  should  consider  this  es¬ 
timate  of  available  suitable  habitat  near  extant  L.  subcoriacea  populations 
and  throughout  the  species'  modeled  distribution  during  its  listing  review. 
However,  further  examination  of  the  strength  of  existing  local  protections, 
the  ability  to  carry  out  management  activities,  and  potential  future  threats 
to  this  suitable  habitat  are  all  necessary. 

Survey  efforts  informed  by  this  SDM  can  be  used  to  improve  sampling  effi¬ 
ciency  and  population  discovery  rates  for  L.  subcoriacea.  However,  imple¬ 
mentation  of  field  surveys  for  L.  subcoriacea  based  on  these  results  will 
likely  benefit  from  knowledge  about  distance  from  roads,  as  well  as  prop¬ 
erty  ownership,  as  extensive  areas  within  multiple  jurisdictions  were  iden¬ 
tified.  Estimates  of  areas  warranting  additional  survey  effort  on  installa¬ 
tions  and  within  their  surrounding  landscapes  can  also  be  used  to  estimate 
costs  for  planning  targeted  searches. 

The  numerous  management  needs  and  opportunities  identified  here  should 
be  pursued  to  the  extent  feasible.  Conservation  partnerships  may  be  helpful 
in  magnifying  the  beneficial  impact  of  these  efforts.  Implementation  will 
likely  reguire  additional  localized  analyses  based  on  criteria  identified  by 
conservation  practitioners  (e.g.,  property  costs  and  acguisition  budgets). 
Additionally,  one  could  consider  other  variables  likely  to  be  important  to 
long  term  L.  subcoriacea  population  success,  including  the  number  of  indi¬ 
viduals  in  populations,  metapopulation  criteria  (e.g.,  size  and  distance  of 
neighbor  populations),  and  ability  to  implement  such  management  actions 
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as  prescribed  fires.  Consequently,  these  results  should  be  viewed  as  an  ini¬ 
tial  screening  to  identify  potential  management  options. 

The  results  of  this  work  indicate  that  the  prospects  for  augmenting  or 
(re)introdudng  L.  subcoriacea  populations  are  good.  However,  before  im¬ 
plementing  any  (re) introductions  it  will  be  important  to  give  consideration 
to  additional  factors  that  can  influence  population  success,  such  as  the  size 
of  the  specific  patches  of  suitable  habitat  and  their  distance  from  neigh¬ 
boring  patches.  Field  validation  of  site  suitability  and  ability  to  implement 
management  actions  are  also  critical  (Guerrant  and  Kaye  2007,  Godefroid 
etal.  2011). 

4.3  Military  Benefits 

Installation  training  ranges  are  essential  for  preparing  DoD  forces  for 
combat  and  complex  missions  across  the  globe.  To  ensure  the  long-term 
sustainability  of  its  training  ranges,  DoD  recognizes  that  it  is  necessary  to 
"sustain  excellence  in  environmental  stewardship"  and  "mitigate  en¬ 
croachment  pressures  on  training  activities"  (USDOD  2015).  Threatened 
and  endangered  species  are  consistently  reported  as  the  primary  encroach¬ 
ment  threat  to  range  accessibility  and  capability  across  DoD,  while  at-risk 
species  are  identified  as  a  key  evolving  challenge  to  sustainability.  The  re¬ 
sults  of  this  effort  support  DoD's  environmental  stewardship  and  en¬ 
croachment  mitigation  requirements  in  multiple  ways. 

Stewardship  of  listed  and  at-risk  species  is  dependent  on  knowing  about 
species'  biology,  locations  where  they  occur,  the  distribution  of  suitable 
habitat,  and  the  viability  of  different  conservation  strategies.  This  effort 
generated  information  related  to  the  latter  three  for  DoD  installations,  re¬ 
gions,  and  across  the  range  of  L.  subcoriacea.  For  DoD  installations,  this 
information  will  improve  survey  efficiency,  aid  development  of  manage¬ 
ment  plans,  expand  available  management  strategies,  and  help  identify 
conservation  partnering  opportunities. 

The  numerous  widely  distributed  L.  subcoriacea  populations  found  across 
installation  landscapes  would  not  only  present  a  management  burden,  but 
also  a  significant  encroachment  constraint  on  DoD  land  use  if  the  species 
were  to  be  listed.  Options  for  mitigating  listed  and  at-risk  species  en¬ 
croachment  on  training  activities  are  dependent  on  species'  listing  status, 
knowledge  about  the  availability  and  distribution  of  suitable  habitat,  and 
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the  viability  of  different  conservation  strategies.  This  effort  generated  in¬ 
formation  related  to  all  of  these  needs.  Increased  efficiency  and  efficacy  of 
new  population  discovery  will  improve  the  accuracy  of  range- wide  esti¬ 
mates  of  L.  subcoriacea  distribution,  population  numbers,  and  total  num¬ 
ber  of  individuals.  Novel  range- wide  estimates  of  available  suitable  habitat 
will  help  identify  the  degree  to  which  L.  subcoriacea  habitats  are  threat¬ 
ened  by  damage  or  destruction.  Both  of  these  are  critical  for  determining 
the  risk  of  extinction  and  listing  status  under  ESA.  This  effort  also  gener¬ 
ated  information  about  landscape-level  conservation  partnering  opportu¬ 
nities  and  areas  potentially  suitable  for  population  (re)introduction,  which 
are  invaluable  strategies  that  DoD  installations  might  wish  to  employ  to 
mitigate  encroachment  on  training  activities. 

The  results  of  this  effort  provide  DoD  natural  resource  managers  with  im¬ 
portant  information  needed  to  adopt  proactive  approaches  to  L.  subcoria¬ 
cea  management,  which  generally  provide  the  greatest  payoff  on  invest¬ 
ment.  Collectively,  the  stewardship  and  encroachment  mitigation  benefits 
generated  will  contribute  towards  improving  the  species'  conservation  sta¬ 
tus,  and  will  help  ensure  that  DoD  maintains  the  flexibility  needed  for  un¬ 
encumbered,  mission- directed  land  use. 
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Appendix  A:  Climate,  Topographic,  Soil,  Fire 
and  Vegetation  Variables  Evaluated 

Of  the  52  climate,  topographic,  soil,  fire,  and  vegetation  variables  evalu¬ 
ated  during  development  of  the  Lindera  suboriacea  habitat  suitability 
map,  only  20  variables  contributing  >  1%  and  having  a  Pearson's  r  <|  0.9| 
were  retained  for  use  in  the  final  model.  Sources  for  the  environmental 
predictors  are:  climate  (www.woridciim.orgl.  fire  regime,  topographic,  and  vege¬ 
tation  (www.landfire.govl.  soils  (www.nrcs.usda.govl. 


Table  A-l.  Fifty-two  climate,  topographic,  soil,  fire  and  vegetation  variables  evaluated  during 
development  of  the  Lindera  suboriacea  habitat  suitability  map. 


Type 

Code 

Name 

Description 

% 

Climate 

Bid 

bioclimatic  indicator  1 

annual  mean  temperature  (°C) 

NA 

Climate 

BI010 

bioclimatic  indicator  10 

mean  temperature  of  warmest  quarter  (°C) 

NA 

Climate 

BI011 

bioclimatic  indicator  11 

mean  temperature  of  coldest  quarter  (°C) 

NA 

Climate 

BI012 

bioclimatic  indicator  12 

annual  precipitation  (mm) 

0.15 

Climate 

BI013 

bioclimatic  indicator  13 

precipitation  of  wettest  month  (mm) 

NA 

Climate 

BI014 

bioclimatic  indicator  14 

precipitation  of  driest  month  (mm) 

NA 

Climate 

BI015 

bioclimatic  indicator  15 

precipitation  seasonality 

3.64 

Climate 

BI016 

bioclimatic  indicator  16 

precipitation  of  wettest  quarter  (mm) 

NA 

Climate 

BI017 

bioclimatic  indicator  17 

precipitation  of  driest  quarter  (mm) 

6.84 

Climate 

BI018 

bioclimatic  indicator  18 

precipitation  of  warmest  quarter  (mm) 

NA 

Climate 

BI019 

bioclimatic  indicator  19 

precipitation  of  coldest  quarter  (mm) 

NA 

Climate 

BI02 

bioclimatic  indicator  2 

mean  monthly  diurnal  temperature  range  (°C) 

NA 

Climate 

BI03 

bioclimatic  indicator  3 

Isothermality  [(BI02/BI07)  *  100]  (%) 

0.82 

Climate 

BI04 

bioclimatic  indicator  4 

temperature  seasonality 

NA 

Climate 

BI05 

bioclimatic  indicator  5 

maximum  temperature  of  warmest  month  (°C) 

NA 

Climate 

BI06 

bioclimatic  indicator  6 

minimum  temperature  of  coldest  month  (°C) 

NA 

Climate 

BI07 

bioclimatic  indicator  7 

temperature  annual  range  [BI05  -  BI06]  (°C) 

NA 

Climate 

BI08 

bioclimatic  indicator  8 

mean  temperature  of  wettest  quarter  (°C) 

NA 

Climate 

BI09 

bioclimatic  indicator  9 

mean  temperature  of  driest  quarter  (°C) 

2.51 

fire  regime 

FRG 

fire  regime  group 

classification  of  presumed  historical  fire  regime 

7.32 

fire  regime 

FRI 

mean  fire  return  interval 

mean  fire  return  interval 

NA 

fire  regime 

PLS 

percentage  of  low 
severity  fire 

percentage  of  low  severity  fire  (%) 

NA 

fire  regime 

PMS 

percentage  of  mixed 
severity  fires 

percentage  of  mixed  severity  fires  (%) 

NA 

fire  regime 

PRS 

percentage  of 
replacement  severity 
fires 

percentage  of  replacement  severity  fires  (%) 

1.77 

fire  regime 


VCC 


vegetation  condition 
class 


classification  of  vegetation  departure  (VDEP) 


12.53 
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Type 

Code 

Name 

Description 

% 

fire  regime 

VDEP 

vegetation  departure 

percentage  estimation  of  vegetation  departure  from 
historical  reference  conditions  (%) 

1.70 

Soils 

AWC 

available  water  capacity 

amount  of  water  than  an  increment  of  soil  depth  can 
store  that  is  available  to  plants 

NA 

Soils 

AWS0150WTA 

available  water  storage 

available  water  storage  to  a  depth  of  150  cm 

NA 

Soils 

CEC7 

cation  exchange  capacity 

amount  of  readily  exchangeable  cations  that  can  be 

NA 

Soils 

CLAY 

percentage  clay 

Soils 

DRAINAGECL 

drainage  class 

Soils 

HYDRICON 

hydric  condition 

Soils 

HYDRICRATI 

hydric  soil 

Soils 

KSAT 

soil  permeability 

Soils 

OM 

organic  matter 

Soils 

PH01MCACL2 

hydrogen  ion  activity  of 
soil 

PH1T01H20 


PTOTAL 

SAND 

SILT 

WTDEPANNMIN 


hydrogen  ion  activity  of 
soil 


phosphorous 
percentage  sand 
percentage  silt 
depth  to  wet  soil  layer 


topographic 

topographic 


topographic 


Vegetation 


Vegetation 

Vegetation 

Vegetation 

Vegetation 

Vegetation 

Vegetation 


WTENTHBAR  field  capacity 

ASPECT  Aspect 

SLOPE  slope  of  terrain 


TPI  topographic  position 

index 


BPS  biophysical  setting 


CC  canopy  cover 

CH  canopy  height 

ESP  environmental  site 

potential 

EVC  existing  vegetative  cover 

EVH  existing  vegetative  height 


electrically  adsorbed  to  negative  charges  in  the  soil, 
soil  constituent,  or  other  material,  at  pH  7.0,  as 
estimated  by  the  ammonium  acetate  method. 


total  clay  (%) 

natural  drainage  conditions  of  the  soil,  refers  to 
frequency  and  duration  of  wet  periods 

natural  condition  of  the  soil  component 
classification  as  hydric  soil  or  not 

amount  of  water  that  would  move  vertically  through  a 
unit  area  of  saturated  soil 

amount  by  weight  of  decomposed  plant  and  animal 
residue 

The  negative  logarithm  to  base  of  10  orthe  hydrogen 
ion  activity  in  the  soil,  using  the  0.01M  CaCI2  method, 
in  a  1:2  soiksolution  ratio.  A  numerical  expression  of 
the  relative  acidity  or  alkalinity  of  a  soil  sample. 

The  negative  logarithm  to  the  base  10,  of  the 
hydrogen  ion  activity  in  the  soil  using  the  1:1  soil- 
water  ratio  method.  A  numerical  expression  of  the 
relative  acidity  or  alkalinity  of  a  soil  sample. 

estimate  of  total  phosphorous  content  of  soil 
total  sand  (%) 
total  silt  (%) 

shallowest  depth  to  a  wet  soil  layer  (water  table)  at 
anytime  during  the  year 

volumetric  content  of  a  soil  water  retained  at  a 
tension  of  1/10  bars  (field  capacity) 

aspect  at  30  meter  resolution 

represents  the  percentage  change  in  elevation 


classified  topographic  position  (e.g.  hilltop,  slope, 
valley) 


vegetation  that  may  have  been  dominant  before  Euro- 
American  settlement 


canopy  cover  (%) 
canopy  height  (m) 

vegetation  that  could  be  potentially  supported  at  a 
site 

vegetation  currently  at  the  site 

mean,  species-over  weighted  existing  vegetation 
height 


existing  vegetation  type  plant  community  classification 
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Appendix  B:  Area  of  L.  subcoriacea  Habitat 
Summarized  by  State,  Ecoregion,  and 
Ownership  for  Different  Management 
Strategies 


Table  B-l.  Total  (km2)  and  percent  (in  parentheses)  area  of  highly  suitable  habitat,  locations 
warranting  additional  survey  effort,  and  locations  potentially  suitable  for  L.  subcoriacea 
population  (re)introductions  across  seven  states  and  associated  ecoregions. 


State,  Ecoregion,  Ownership 

Total 

Suitable 

Additional 

Survey 

Reintroduction 

Alabama 

94219.4  (14.7) 

1870.7  (16.1) 

1118.7  (23.7) 

250.2  (27.6) 

Coastal  Gulf  Plains 

35878.5  (5.6) 

432.6(3.7) 

157.2  (3.3) 

5.3  (0.6) 

Federal 

498.9  (0.1) 

1.4  (0) 

1.3  (0) 

1.3  (0.1) 

NGO 

6.2  (0) 

0.3  (0) 

0.3  (0) 

0.3  (0) 

Other 

582.6  (0.1) 

16.3  (0.1) 

0.9  (0) 

0.9  (0.1) 

Private 

34129.8  (5.3) 

407.1  (3.5) 

151.4  (3.2) 

0(0) 

State 

661  (0.1) 

7.6  (0.1) 

3.3  (0.1) 

2.9  (0.3) 

Gulf  Plains 

45630.6  (7.1) 

1363.2  (11.7) 

928.3  (19.7) 

240.1  (26.5) 

Federal 

713.6  (0.1) 

229.4  (2) 

221.7  (4.7) 

220.8  (24.3) 

NGO 

1.9  (0) 

0(0) 

0(0) 

0(0) 

Other 

591.9  (0.1) 

17.9  (0.2) 

14.4  (0.3) 

14.1  (1.6) 

Private 

43745.3  (6.8) 

1109.6  (9.5) 

686.7  (14.6) 

0(0) 

State 

577.9  (0.1) 

6.4  (0.1) 

5.5  (0.1) 

5.2  (0.6) 

Southeastern  Coastal  Plain 

226.3  (0) 

0.4  (0) 

0.1  (0) 

0.1  (0) 

NGO 

0(0) 

0(0) 

0(0) 

0(0) 

Other 

93.9  (0) 

0.1  (0) 

0.1  (0) 

0.1  (0) 

Private 

131.9  (0) 

0.2  (0) 

0(0) 

0(0) 

State 

0.5  (0) 

0(0) 

0(0) 

0(0) 

Southern  Piedmont 

12484.1  (2) 

74.6  (0.6) 

33.1  (0.7) 

4.6  (0.5) 

Federal 

776.5  (0.1) 

0.4  (0) 

0.3  (0) 

0.3  (0) 

NGO 

1.9  (0) 

0(0) 

0(0) 

0(0) 

Other 

201  (0) 

4.5  (0) 

2.3  (0) 

2.3  (0.2) 

Private 

11320  (1.8) 

65.6  (0.6) 

28.3  (0.6) 

0(0) 

State 

184.7  (0) 

4.2  (0) 

2.2  (0) 

2.1  (0.2) 

Florida 

71575.4  (11.2) 

720.6(6.2) 

147.3  (3.1) 

56.1(6.2) 

Coastal  Gulf  Plains 

20600.4  (3.2) 

710.3  (6.1) 

145  (3.1) 

55.8(6.2) 

Federal 

129.4  (0) 

0(0) 

0(0) 

0(0) 

NGO 

26.7  (0) 

1.3  (0) 

1.2  (0) 

1.2  (0.1) 

Other 

3082.3  (0.5) 

56.7  (0.5) 

37.6  (0.8) 

37.3  (4.1) 
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State,  Ecoregion,  Ownership 


Private 

State 

Southeastern  Coastal  Plain 


Total 


16160.3  (2.5) 

1201.8  (0.2) 
50975  (8) 


Georgia 
Gulf  Plains 
Other 
Private 

Southeastern  Coastal  Plain 


137737.7  (21.5) 
3.9  (0) 
0.8  (0) 
3.1(0) 

93155  (14.6) 


Federal 

2097.7  (0.3) 

NGO 

636  (0.1) 

Other 

4019.1  (0.6) 

Private 

84750  (13.3) 

State 

1652.2  (0.3) 

Southern  Piedmont 
Federal 
NGO 
Other 
Private 
State 

Louisiana 

Coastal  Gulf  Plains 
Federal 
NGO 
Other 
Private 
State 

Mississippi  Delta 


44578.9  (7) 

834.7  (0.1) 

52.2  (0) 

1920  (0.3) 

41277.6  (6.5) 

494.3  (0.1) 

53731.5  (8.4) 
11298  (1.8) 

145.8  (0) 

16.8  (0) 

82.6  (0) 

10796.1(1.7) 

256.7  (0) 

42433.4  (6.6) 


Suitable 


629.7  (5.4) 

22.6  (0.2) 

10.3  (0.1) 


Federal 

4896.6  (0.8) 

0(0) 

NGO 

89.8  (0) 

0.1(0) 

Other 

4615  (0.7) 

0.2  (0) 

Private 

36928  (5.8) 

9.3  (0.1) 

State 

4445.5  (0.7) 

0.7  (0) 

238.5  (2.1) 
0(0) 
0(0) 
0(0) 

151.6  (1.3) 
0.1  (0) 
0.9  (0) 

6.4  (0.1) 

140.7  (1.2) 

3.5  (0) 

86.9  (0.7) 

6.5  (0.1) 
0(0) 

13.6  (0.1) 

61.4  (0.5) 

5.4  (0) 

1794.6  (15.4) 

1657.4  (14.3) 

1.8  (0) 
0.1(0) 

34.2  (0.3) 
1588.8(13.7) 

32.4  (0.3) 

137.2  (1.2) 


Additional 

Survey 


88.8(1.9) 

17.5  (0.4) 
2.2  (0) 
0(0) 
0(0) 
0.1  (0) 

1.9  (0) 
0.2  (0) 

99.5  (2.1) 
0(0) 
0(0) 
0(0) 

70.8(1.5) 
0.1  (0) 
0.4  (0) 

3.2  (0.1) 
64  (1.4) 

3.1  (0.1) 

28.7  (0.6) 

4.7  (0.1) 

0(0) 

7.7  (0.2) 

12.4  (0.3) 

3.9  (0.1) 

634.6  (13.4) 

532.8  (11.3) 
1.8  (0) 
0.1  (0) 
1.6  (0) 

507.1  (10.7) 

22.2  (0.5) 

101.7  (2.2) 


Reintroduction 


0(0) 

17.4  (1.9) 
0.3  (0) 
0(0) 
0(0) 
0.1  (0) 
0(0) 
0.2  (0) 

23  (2.5) 
0(0) 
0(0) 
0(0) 

6.7  (0.7) 
0.1  (0) 
0.4  (0) 

3.2  (0.4) 
0(0) 

3.1  (0.3) 

16.2  (1.8) 

4.7  (0.5) 

0(0) 

7.6  (0.8) 
0(0) 

3.9  (0.4) 

37.8  (4.2) 

25.3  (2.8) 

1.8  (0.2) 
0.1  (0) 

1.5  (0.2) 
0(0) 

21.8  (2.4) 

12.5  (1.4) 


Federal 

792.2  (0.1) 

0.1(0) 

0.1  (0) 

0.1  (0) 

NGO 

0.7  (0) 

0(0) 

0(0) 

0(0) 

Other 

610.3  (0.1) 

1.3  (0) 

1.2  (0) 

1.2  (0.1) 

Private 

38027.5  (5.9) 

122.5  (1.1) 

89.1(1.9) 

0(0) 

State 

3002.8  (0.5) 

13.3  (0.1) 

11.3  (0.2) 

11.2  (1.2) 
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State,  Ecoregion,  Ownership 


Mississippi 
Coastal  Gulf  Plains 
Federal 


Eastern  Coastal  Plain 
Federal 
NGO 
Other 
Private 
State 

Northeastern  Piedmont 
Federal 
NGO 
Other 
Private 


Total 


103844  (16.2) 
33416.6  (5.2) 

2597.9  (0.4) 
50  (0) 


Suitable 


1598.8(13.7) 

1042.4  (9) 

138.7  (1.2) 

7.4  (0.1) 


Additional 

Survey 


1076.9  (22.8) 

767.9  (16.3) 

132.8  (2.8) 

7.1  (0.1) 


51528.2  (8.1) 

2134.7  (0.3) 

57.4  (0) 

1822.7  (0.3) 
44259  (6.9) 

3254.4  (0.5) 

46049.1  (7.2) 
226.5  (0) 

50.7  (0) 

695.4  (0.1) 

44569.4  (7) 


120.4  (1) 

2.2  (0) 
0(0) 
1-7(0) 

108.4  (0.9) 

7.9  (0.1) 

4850.2  (41.7) 

76.8  (0.7) 

8.8  (0.1) 

48.7  (0.4) 

4616.3  (39.7) 


39.2  (0.8) 

1.1  (0) 
0(0) 
0.5  (0) 

31.8  (0.7) 

5.8  (0.1) 

1482.2  (31.4) 

73.2  (1.6) 

6.4  (0.1) 

9.4  (0.2) 
1319  (27.9) 


Reintroduction 


344.7  (38) 

287.9  (31.8) 

131.9  (14.5) 
7  (0.8) 


Other 

858.6  (0.1) 

29.3  (0.3) 

25.1  (0.5) 

24  (2.6) 

Private 

28479.5  (4.5) 

732  (6.3) 

477.4  (10.1) 

0(0) 

State 

1430.7  (0.2) 

134.9  (1.2) 

125.6  (2.7) 

125  (13.8) 

Gulf  Plains 

69928.3  (10.9) 

556.2  (4.8) 

308.9  (6.5) 

56.7  (6.3) 

Federal 

2174.5  (0.3) 

42.1  (0.4) 

35.8(0.8) 

35.5  (3.9) 

NGO 

0.3  (0) 

0(0) 

0(0) 

0(0) 

Other 

936.5  (0.1) 

17.3  (0.1) 

13.5  (0.3) 

13.4  (1.5) 

Private 

65331.6  (10.2) 

482.3  (4.1) 

251.8  (5.3) 

0(0) 

State 

1485.3  (0.2) 

14.5  (0.1) 

7.9  (0.2) 

7.8  (0.9) 

Mississippi  Delta 

499.1  (0.1) 

0.2  (0) 

0.1  (0) 

0(0) 

Federal 

69  (0) 

0(0) 

0(0) 

0(0) 

Other 

78.7  (0) 

0(0) 

0(0) 

0(0) 

Private 

336.8  (0.1) 

0.1(0) 

0.1  (0) 

0(0) 

State 

14.6  (0) 

0(0) 

0(0) 

0(0) 

North  Carolina 

99034.6  (15.5) 

4977.1  (42.8) 

1522.4  (32.3) 

167.5  (18.5) 

6.6  (0.7) 
1.1  (0.1) 
0(0) 
0.4  (0) 
0(0) 

5.1  (0.6) 
160.5  (17.7) 

72.5  (8) 

6.1  (0.7) 
8.9  (1) 
0(0) 


State 

507.2  (0.1) 

99.6  (0.9) 

74.2  (1.6) 

73.1(8.1) 

Southern  Piedmont 

1457.2  (0.2) 

6.5  (0.1) 

1.1  (0) 

0.4  (0) 

NGO 

1.1  (0) 

0(0) 

0(0) 

0(0) 

Other 

14  (0) 

0.1  (0) 

0(0) 

0(0) 

Private 

1432.7  (0.2) 

6  (0.1) 

0.7  (0) 

0(0) 

State 

9.4  (0) 

0.4  (0) 

0.4  (0) 

0.4  (0) 
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State,  Ecoregion,  Ownership 


South  Carolina 

Eastern  Coastal  Plain 
Federal 
NGO 
Other 
Private 
State 

Northeastern  Piedmont 
Federal 
NGO 
Other 
Private 
State 

Southeastern  Coastal  Plain 
Federal 
NGO 
Other 
Private 
State 

Southern  Piedmont 
Federal 
NGO 
Other 
Private 
State 

Grand  Total 


Total 


79402.4  (12.4) 

50635.2  (7.9) 

1555.6  (0.2) 

64.9  (0) 

2717.6  (0.4) 

44788.2  (7) 

1508.8  (0.2) 

4696.7  (0.7) 

16.6  (0) 

0(0)  ~ 

20.5  (0) 

4605.5  (0.7) 

54.1  (0) 

847.7  (0.1) 

63.5  (0) 

0(0)  ~ 

144  (0) 

543.2  (0.1) 

97  (0) 

23222.8(3.6) 

1191.8  (0.2) 

2.9  (0) 

595.2  (0.1) 

21245.8  (3.3) 
187  (0) 

639545.1  (100) 


Suitable 


428.3  (3.7) 

65.5  (0.6) 

0.8  (0) 

0(0) 

4.8  (0) 

57.7  (0.5) 

2.2  (0) 

288.9  (2.5) 

5.5  (0) 

0(0) 

4.4  (0) 

264.7  (2.3) 

14.3  (0.1) 

0.1  (0) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

73.9  (0.6) 
6.6  (0.1) 

0(0) 

2.9  (0) 

63.5  (0.5) 

1(0) 

11628.6  (100) 


Additional 

Survey 


119.9  (2.5) 
23.1(0.5) 

0-7  (0) 

0(0) 

3  (0.1) 

17.7  (0.4) 

1.6  (0) 

85  (1.8) 

5.5  (0.1) 
0(0) 

0-7  (0) 

65.7  (1.4) 
13.1(0.3) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

11.8  (0.2) 

1.4  (0) 

0(0) 

1.2  (0) 

8.5  (0.2) 

0-7  (0) 

4719.3  (100) 


Reintroduction 


27.6  (3) 

5.1  (0.6) 

0.7  (0.1) 

0(0) 

2.7  (0.3) 
0(0) 

1.6  (0.2) 

19.3  (2.1) 

5.5  (0.6) 

0(0) 

0.7  (0.1) 
0(0) 

13.1  (1.4) 
0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

3.3  (0.4) 

1.4  (0.2) 
0(0) 

1.2  (0.1) 
0(0) 

0.7  (0.1) 

906.9  (100) 
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Acronyms  and  Abbreviations 


Term 

Definition 

AFB 

Air  Force  Base 

ASA(ALT) 

Assistant  Secretary  of  the  Army  for  Acquisition,  Logistics,  and 
Technology 

AUC 

Under  the  Curve 

BPS 

Biophysical  Settings 

CART 

Classification  and  Regression  Tree 

CC 

Canopy  Cover 

CEERD 

U.S.  Army  Corps  of  Engineers,  Engineer  Research  and  Development 
Center 

CERL 

Construction  Engineering  Research  Laboratory 

CH 

Canopy  Height 

DoD 

U.S.  Department  of  Defense 

DOE 

U.S.  Department  of  Energy 

DOI 

U.S.  Department  of  Interior 

ERDC 

U.S.  Army  Engineer  Research  and  Development  Center 

ERDC-CERL 

Engineer  Research  and  Development  Center,  Construction  Engineering 
Research  Laboratory 

ESA 

U.S.  Endangered  Species  Act 

ESP 

Environmental  Site  Potential 

ESRI 

Environmental  Systems  Research  Institute,  Inc. 

FRG 

Fire  Regime  Group 

GAP 

USGS  Gap  Analysis  Program 

JFTC 

Joint  Forces  Training  Center 

LANDFIRE 

Landscape  Fire  and  Resource  Management  Planning  Tools 

NGO 

Nongovernmental  Organization 

NLCD 

National  Land  Cover  Data 

NRCS 

Natural  Resources  Conservation  Service 

OSI 

Open  Systems  Interconnection 

PRS 

Percentage  of  Replacement  Severity  (fires) 

ROC 

Receiver  Operating  Curve 

SAR 

Species  at  Risk 

SDM 

Species  Distribution  Modeling 

SF 

Standard  Form 

SSURGO 

Soil  Survey  Geographic 

TM 

Thematic  Mapper 

TNC 

The  Nature  Conservancy 

TPI 

Topographic  Position  Index 

TR 

Technical  Report 

URL 

Universal  Resource  Locator 
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Definition 

U.S. 

United  States 

USA 

United  States  of  America 

USFWS 

U.S.  Fish  and  Wildlife  Service 

USGS 

U.S.  Geological  Survey 

VCC 

Vegetation  Condition  Class 

VDEP 

Vegetation  Departure 

WTDEPANMIN 

Lower  Depth  to  the  Water  Table 

WTDEPANNMIN 

Depth  to  Wet  Soil  Layer 

WTENTHBAR 

Field  Capacity 
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